Electron spin decoherence in quantum dots due to interaction with nuclei 
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We study the decoherence of a single electron spin in an isolated quantum dot induced by hyperfine 
interaction with nuclei for times smaller than the nuclear spin relaxation time. The decay is caused 
by the spatial variation of the electron envelope wave function within the dot, leading to a non- 
uniform hyperfine coupling. We evaluate the spin correlation function with and without magnetic 
fields and find that the decay of the spin precession amplitude is not exponential but rather power 
(inverse logarithm) law-like. For fully polarized nuclei we find an exact solution and show that the 
precession amplitude and the decay behavior can be tuned by the magnetic field. The corresponding 
decay time is given by hN/A, where yl is a hyperfine interaction constant and A'^ the number of nuclei 
inside the dot. The amplitude of precession, reached as a result of the decay, is finite. We show that 
there is a striking difference between the decoherence time for a single dot and the dephasing time 
for an ensemble of dots. 
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The spin dynamics of electrons in semiconducting 
nanostructures has become of central interest in recent 
years Q]. The controlled manipulation of spin, and in 
particular of its phase, is the primary prerequisite needed 
for novel applications in conventional computer hardware 
as well as in quantum information processing. It is thus 
desirable to understand the mechanisms which limit the 
spin phase coherence of electrons, in particular in GaAs 
semiconductors, which have been shown to exhibit un- 
usually long spin decoherence times T2 exceeding 100 ns 
1^ . Since in GaAs each nucleus carries spin, the hyperfine 
interaction between electron and nuclear spins is unavoid- 
able, and it is therefore important to understand its effect 
on the electron spin dynamics. This is particularly so for 
electrons which are confined to a closed system such as 
a quantum dot with a spin 1/2 ground state, since, be- 
sides fundamental interest, these systems are promising 
candidates for scalable spin qubits [^. For recent work 
on spin relaxation (characterized by Ti times) in GaAs 
nanostructures we refer to Refs. 



Motivated by this we investigate in the following the 
spin dynamics of a single electron confined to a quantum 
dot in the presence of nuclear spins. We treat the case of 
unpolarized nuclei perturbatively, while for the fully po- 
larized case we present an exact solution for the spin dy- 
namics and show that the decay is non-exponential and 
can be strongly influenced by external magnetic fields. 
We use the term " decoherence" to describe the case with 
a single dot, and the term "dephasing" for an ensemble 
of dots 1^ . The typical fluctuating nuclear magnetic field 
seen by the electron spin via the hyperfine interaction is 
of the order of Q '--^ A/ {^/Ngfis), with an associated 
electron precession frequency lon c± A/ ^/N, where A is 
a hyperfine constant, g the electron g-factor, and fiB the 
Bohr magneton. For a typical dot size the electron wave 
function covers approximately N = 10^ nuclei, then this 



field is of the order of 100 Gauss in a GaAs quantum dot. 
The nuclei in turn interact with each other via dipolar 
interaction, which does not conserve the total nuclear 
spin and thus leads to a change of a given nuclear spin 
configuration within the time T„2 ~ 10~*s, which is just 
the period of precession of a nuclear spin in the local 
magnetic field generated by its neighbours. 

We note that there are two different regimes of inter- 
est, depending on the parameter lonTc where Tc is the 
correlation time of the nuclear field. The simplest case 
is given by the perturbative regime ujnTc <C 1, charac- 
terized by dynamical narrowing: different random nu- 
clear configurations change quickly in time, and, as a 
result, the spin dynamics is diffusive with a dephasing 
time ~ l/(w^rc). A more difficult situation arises when 
^nTc 3> 1, requiring a nonperturbative approach. It is 
this regime which we will consider in this paper, i.e. the 
electron is localized in a quantum dot, and the correla- 
tion time is due to the internal nuclear spin dynamics, 
i.e., Tc = T„2, giving ojmTc — lO'*. Next, we need to ad- 
dress the important issue of averaging over different nu- 
clear spin configurations in a single dot. Without internal 
nuclear spin dynamics, i.e. T„2 — > 00, no averaging is in- 
dicated. However, each flip-flop process (due to hyperflne 
interaction) creates a different nuclear configuration, and 
because of the spatial variation of the hyperfine coupling 
constants inside the dot, this leads to a different value 
of the nuclear field seen by the electron spin and thus to 
its decoherence. Below we will find that this decoherence 
is non-exponential, but still we can indicate a character- 
istic time given by {A/hN)~^ 0. Moreover, we shall 
find that T„2 > {A/hN)~^, and thus stiU no averaging 
over the nuclear configurations is indicated (and dipolar 
interactions will be neglected henceforth). To underline 
the importance of this point, we will contrast below the 
unaveraged correlator with its average. 
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Unpolarized nuclei. We consider a single electron con- 
fined to a quantum dot whose spin S couples to an exter- 
nal magnetic field B and to nuclear spins {!'} via hyper- 
fine contact interaction, described by the Hamiltonian 



7i = SMsS • B + S • h 



N ■ 



(1) 



where h.N/{g^B) is the nuclear field. The coupling con- 
stant with the i-th nucleus, Ai = Auo|^'(ri)p, contains 
the electron envelope wave function ^'(ri) at the nuclear 
site Ti, and vq is the volume of the crystal cell. We start 
with the case B = 0, and for simplicity we consider nu- 
clear spin 1/2. Neglecting dipolar interactions between 
the nuclei, we can consider only some particular nuclear 
configuration, described in the eigenbasis as \{Il} >, 
with II = ±1/2. Moreover, we assume an unpolarized 
configuration with a typical net nuclear magnetic field 
A/ {\/N g^B), being much less than A/ (fully polar- 
ized case) . We study now the decay of the electron spin 
from its initial [t = 0) 5z-eigenstate | 'f|->. For this we 
evaluate the correlator C„(t) n\5Sz{t)Sz\n >, where 
5Sz{t) = Sz{t) - Sz, and S'^(i) exp(it7T!)S'^ exp(-ii7i:). 
Since at < = the total (electron and nuclear) state 
W >— I tl", > is an eigenstate of Tio = SzHnz 
(with eigenergy e„), we can expand in the perturbation 

V = {l/2){S+hN-+S-hN+) {withn = no + V)- Going 
over to the interaction picture, we obtain in leading order 



Cn{t) 



\Vnk 



k 



-{C0s{ujnkt) - 1) 



(2) 



where Vnk =< n\V\k > is the matrix element between 
initial state n —'^,{...,1^ = —1/2,...} and intermediate 
state k =-11, {...., = +1/2, ...}, and w„fc = e„ — Cfe. Us- 
ing iKfeP = Al < n\l/2 - i^\n > /4, and Lo.^k = (/i.)„ + 
Ak/2, where {hz)n =< n\hNz\n >, we obtain for the typ- 
ical nuclear configuration, for which {hz)n — '^n ^ 



C„(<)~ 7 [lo - hir) cos{{hz)nt) + hir) sm{{hz)nt)] 



dzxUz)F,iTxUz)), 0,1,2, 



(3) 



where 7 = -Ay{8nN{hz)l), Fo = 1/2, ^1(77) = 
sin 77/7/ + (cos?7 — l)/?7^, and F2{ri) = sin 77/77^ — cosrj/rj. 
Here, N — azO? /vq 3> 1 is the number of nu- 
clei inside the dot, and r = At/2TTN. We have 
replaced the sums over k (which run over the en- 
tire space) by integrals and used that |^(p, z)^ — 
(l/Tra^a^) exp(— p^/a^)Xo(^)- Here a,az are the dot 
sizes in the lateral and transverse (perpendicular to 
the 2D plane) directions, resp., and the transverse 
wave function is normalized, i.e. dzx^iz) = 1. 

For any analytic function with expansion Xo(-^) ~ 
Xo(0) ~ •2^^(Xo)"/2 near its maximum, we have /i.2(t ^ 
1) = ±(xg(0)/r3/2)y^f7(^[sin(rxg(0))Tcos(rxg(0))]. 



Thus, Cn{t) oc 1/t3/2, i.e. the spin decay follows a 
power law for times r 3> 1, i.e. t :s> {A/N)~^. Note 
that for the typical nuclear configuration the quantity 

/N{hz)n is of order unity, thus the part of the elec- 
tron spin state which decays is of the order of the initial 
value. Hence the same holds for the spin part which sur- 
vives at r ^ 1. For the fully polarized nuclear state the 
result (||) should be multiplied by 2, and (/iz)„ should 
be replaced by A/ 2. Moreover, in the presence of a large 
Zeeman field, = g^J^uBz 3> u>n, we should substitute 
tz for {hz)n- Thus, the asymptotic behaviour of C„(t) is 
not changed, the only difference being that the decaying 
part of the initial spin state is small now, being of the 
order of (wat/ez)^ <C 1. 

We note that C„(t) in (||) is quasiperiodic in t, and, 
thus, it will decay only up to the Poincare recurrence time 
Tp. This time can be found from the condition that the 
terms omitted when converting sums to integrals become 
comparable with the integral itself. This will happen at 
T ~ iV, giving Tp = 0.1s — Is. 

In next order, F^, we face the problem of " resonances" , 
i.e. the corrections will contain zero denominators. This 
gives rise to linearly growing terms cx LONt, even for 
t ^ [A/N)^^. In higher order the degree of the di- 
vergence will increase. This means that the decay law 
we found can, in principle, change after proper resumma- 
tion, because no small expansion parameter exists, which, 
strictly speaking, would justify a perturbative approach. 
Still, the result found in lowest order remains qualita- 
tively correct in that it shows that a non-uniform hy- 
perfine coupling leads to a non-exponential decay of the 
spin. This conclusion is confirmed by an exactly solvable 
case to which we turn next. 

Polarized nuclei. Exactly solvable case. In this sec- 
tion we consider the exactly solvable case where the ini- 
tial nuclear spin configuration is fully polarized. We 
also allow for a magnetic field but neglect its effect 
on the nuclear spins. With the initial wave function 
^0 = I ^sTjTiT:--- > we can construct the exact wave 
f miction of the system for t > {), 

■^{t) = «(*)*„ + Wl T, T, ik, T, ... >, (4) 

k 

with normalization |Q!(t)P +X]fe l/3fe(i)P = Ij and we as- 
sume that a{t = 0+) — l,a{t < 0) = 0. Then, inserting 
^(t) into the Schrodinger equation we obtain 



.da{t) _ 1 
' dt ~ ~4 



A, 



Aa{t)+J2^Mt) 



2 ' 



.^ = (^-^)A(i) + ^«(t) + ^, (5) 
dt M 2 ^ ^ 2 ^ ^ 2 ' ^ ' 

where A ~ Ak- Laplace transforming (^, we obtain 
If dioi exp[(w - iA'/4:)t] 



ait) 



2TTi 



ito + ez + irNiuj J dz ln(l 



2ttNu ' 



(6) 
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where A' = A + 2ez. We have replaced the sum 
Sfe - -^A'/4+A 12 ^^^-"^ xy-plane by an integral and 



i(^-A'/4+Afc/2 

calculated it using |\E'(rfc)p given above. As usual, the 
integration contour F in Eq.(|^) is the vertical line in 
the complex w-plane so that all singularities of the inte- 
grand lie to its left. These singularities are: two branch 
points (lo = O,cjo = *^Xo(0)/27rA^), and first order 
poles which lie on the imaginary axis (lo = iv). For 
Ez > there is one pole, while for ez < there are 
two poles, and for — Q there is one first order pole 
at LOi K, iA/2 + zv4 J dzxt{z) / ^nN . For the contribution 
from the branch cut between w = and lo — loq wc obtain 



-iA't/i 

ttN 



dK2zoKe^^''^^[K j dzln|l 



(7) 



where r' = rxo(O), and zq = zq{k) is defined through 
Xo(-zo) = Xo(O)'*- Let us consider first the case ez = 0. 
The asymptotic behavior of the integral (|^) for r ^ 1 is 
determined by k ^ 1. For example, for Xssi^) / Xoi^) — 
exp(— z^) we find a cx 1/ In^^^ r. Thus, the decay of |a(i)| 
starts at r > 1, i.e. at i > A^^N , as in the unpolarized 
case. Note that the magnitude of a is of order 1/iV, thus 
the decaying part of the initial spin state has this small- 
ness as well, in contrast to the unpolarized case above 
where this part is of order unity For large Zeeman 
field (|ez| ^ A) and for r ^ 1, the main contribution in 
is given for n ~* 1. Expanding XoC-^) ^'^^ small z (see 
above), we obtain Zq — 2x§(0)(l — n)/{Xo)" ■ Then from 
Eq.(0) we have for |ez| > A 



a(T > 1) 



(8) 



From this we find then that the correlator Co(i) = — < 
■^o\5Sz{t)Sz\^Q >= (1 - |a(t)P)/2 agrees with the per- 
turbative result (||) for the fully polarized state, i.e. 
Co(i) c>c Xjt?'!'^. This agreement is to be expected, since 
for large Zeeman field, the perturbative treatment with 
a small parameter Al\ez\ <C 1 is meaningful How- 
ever, at zero Zeeman field, when the system cannot be 
treated perturbatively, we find Co(t) oc 1/ In'^^^ and the 
agreement with (||) breaks down. Nevertheless, the char- 
acteristic time scale for the onset of the non-exponential 
decay is the same for all cases and given by (A/iV)"^. 

There are several interesting features which we can ob- 
serve for the fully polarized state. In an external Zee- 
man field, the effective gap seen by the electron spin 
is A 12 = yl/2 + e^. Thus, when tz is made negative 
this gap decreases and even vanishes at \tz\ = A/2. 
From Eq.(^ we find that the two poles are symmetric 
in this case, and the system resonates between the two 
frequencies w± = ±iA(/ xt{z)dzYl'^ j \/%-nN . Note that 
the residual gap is of order A/^fN (and not A/N, as 



one might naively expect). Near this Zeeman field we 
have |a(t)p = cos^(u;+t) (up to small corrections of or- 
der 1/A^), and, as a result, |ap averaged over time is 
1/2, i.e. the up and down states of the electron spin are 
strongly coupled via the nuclei. In contrast, outside this 
resonance regime the value of |ap is close to 1 (again 
with small 1/A^ corrections), i.e. < Sz{t) >= 1/2 - jap 
is close to -1/2 at any time. The width of the resonance 
is ~ A/\/N, i.e. small compared to the initial gap A/2. 
We note that this behavior represents periodic (Rabi) os- 
cillations with a single well-defined frequency and is not 
related to decoherence. [The latter is described by the 
branch cut contribution d which remains small (order 
1/N) even near the resonance.] This abrupt change in 
the amplitude of oscillations of < Sz{t) > (when chang- 
ing ez in a narrow interval around A/2) can be used for 
an experimental detection of the fully polarized state. 
Note that the weight of the upper pole alone (i.e. that 
which exists at ez = 0) also drops abruptly from a value 
close to 1 to a value much smaller than 1 in the same 
narrow interval, which can be experimentally checked by 
Fourier analysis. Another special value of Zeeman field 
corresponds to the case when the upper pole is close to 
ujQ (k = 1)- the upper edge of the branch cut. This oc- 
curs (see Eq.(^)) at the critical value e* = bA/2 < 0, 
where b = xl{0) J dz In \1 - xl{z)/Xom < -1 is a 
non-universal number which depends on the dot shape. 
Since at finite Zeeman field the asymptotics in t is de- 
termined by k's close to 1, we see from Eq.(^ that for 
e « ej the asymptotics changes abruptly. Indeed, for 
{{ez - e*z)/Af < 1, we find d oc l/Vr, for 1 < t < 
((ez-e*)/A)-2, and d oc l/r^/^, for r » {{ez~e*)/ A)-^ . 
Thus, when approaching the critical Zeeman field e\ 
there is a slow down of the asymptotics from I/t^^^ to 
l/r^/^. It is interesting that this slow down is related 
to a strong modulation of the density of states (DOS) of 
the excitations within the continuum band (branch cut) 
near its edge when — > e*. In the subspacc of none 
or one nuclear spin flipped (see Eq.(^), the DOS be- 
comes ^{u) — Im[Go{u) + d/duhiD{u)], where u — ico, 
Go{u) — l/{u + Ak/2) is the "unperturbed Green's 
function", and D{u) is the denominator of a{Lu), see 
Eq.(||). One can then show that for e* (i.e. the 

upper pole approaches the continuum edge), the DOS 
develops a square root singularity: v{u) oc I/^/wq — u. 
Simultaneously, the weight of the upper pole vanishes 
linearly in ez as e* — ez — > 0. 

Finally, the nuclear state is characterized by Pk{t), 
which allows for similar evaluation as a. Here we just 
note that its branch cut part, Pk{t), is nonmonotonic in 
time, particularly pronounced at ez — > e*: First (3k{t) 
grows like y^r, until r reaches ~ 1/(1 — a^) ^ 1, and 
then it decays like 1/(^/7(1 — Ofe)), with Ofe = Ak/Ao 1. 
Thus, Pk is maximal for Ak close to Aq, i.e. the nuclei 
near the dot-center are affected most by the hyperfine 
interaction with the electron spin. 
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Averaging over nuclear configurations. We have seen 
that the decay of Cn(t) occurs in the time interval 
N/A < t < N'^/A, with N/A ~ IQ-^s in GaAs dots. 
On the other hand, the electron spin precesses in the net 
nuclear field (see Eq. (^) with the characteristic period 
~ Wj^i ~ 10-^ ^ IQ-^s. Thus, < N/A, and 
we see that the electron spin undergoes many precessions 
in a given nuclear field before decoherence sets in due to 
the non- uniform hyperfine couplings Ak- This behavior 
changes dramatically when we average over nuclear con- 
figurations [0. For that purpose we consider high tem- 
peratures, ksT hcjM, and average Cn{t) in Eq.(||) over 
all nuclear configurations, i.e. C{t) = ^nCn{t)/J2n- 
We then find 



C{t) 



For 



dh / dt2Y[ 



r '^'i / 
COS[ (^1 



t2)]. (9) 



T « 1, we get ni5^fcCOs(Aii/2) = 
exp[^NC{At/2TTN)^], where C = tt / dzxt{z)/4:. Thus, 
the averaged spin correlator C{t) is of order 

— J^"* dx^{x), with <I> being the error function. Thus, 
C{t) grows without bound as wjv^ for cjjv^ ^ 1 (the 
condition r <C 1 can still be satisfied). Consequently, the 
perturbative approach breaks down even in leading order 
in V (we recall that without averaging the divergences 
occur in all higher but not in lowest order). To treat 
this case properly, we need a non-perturbative approach. 
For that purpose we calculate now the correlator C{t) 
exactly by treating the nuclear field purely classically, 
i.e. as a c-number. Then we obtain. 



Cnit) 



/)2 

4h% 



(1 — cos hj^t), 



(10) 



where h 



"'N± 



N — 
"'Nx 



5~ 

Nz 

"■Ny 



hjf^ is the nuclear field. 



with 

The value of hff corresponds to 
a given nuclear configuration n. To make contact with 
the perturbation procedure which we used before in the 
quantum case we go to the regime h% ^ <C hj^^^, where 
/lAT can be replaced by h^z in Eq.(|lQ). Then we av- 
erage the resulting expression ihj^j_/hlj^){l — cosh^zt) 
over a Gaussian distribution for h^, i.e. over P{hj\j) cx 
exp(— 3/i^/2a;^). The result becomes proportional to 
dzexp(— z^/2)(l — cos(7z))/z^ oc dx^{x), where 
7 — cjjvVVS- Thus, we see that we obtain exactly the 
same functional form as before from Eq. (|^) with the same 
divergencies in t. This reassures us that the treatment of 
the nuclear field as a classical field is not essential. On 
the other hand, the same Gaussian averaging procedure 



can now be applied to the non-perturbative form Eq.(lO) 
Defining Cci{t) — j dhfqP[ht^)Cn{t)j we obtain 



Ccl{t) 



(11) 



Thus we get rapid (Gaussian) decay of the correlator for 
t ^ w^^, giving the dephasing time = y/N /A. This 



means that < Sz{t)Sz > saturates at 1/3 of its initial 
value of 1/4. Finally, it seems likely that for the case of 
nuclear quantum spins a non-perturbative treatment of 
the averaged correlator C{t) will lead to a similar rapid 
time decay as found for the classical case in Eq. (pi]). 

In conclusion, we have studied the spin decoherence 
of an electron confined to a single quantum dot in the 
presence of hyperfine interaction with nuclear spins. The 
decoherence is due to a non-uniform coupling of the elec- 
tron spin to nuclei located at different sites. The decoher- 
ence time is given by UN/ A and is of the order of several 
IIS. It is shown that in a weak external Zeeman field the 
perturbative treatment of the electron spin decoherence 
is impossible. Moreover, the decay of the electron spin 
correlator in time does not have an exponential charac- 
ter, instead it is given by a power or inverse logarithm 
law. We have shown that there is a strong difference be- 
tween the decoherence time for a single dot, TiN/A, and 
the dephasing time for an ensemble of dots, h^/N /A. 
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